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Abstract. We present a mathematical proof of Einstein's formula for the effective viscosity of 
a dilute suspension of rigid neutrally-buoyant spheres when the spheres are centered on the vertices 
of a cubic lattice. We keep the size of the container finite in the dilute limit and consider boundary 
effects. Einstein's formula is recovered as a first-order asymptotic expansion of the effective viscosity 
in the volume fraction. To rigorously justify this expansion, we obtain an explicit upper and lower 
bound on the effective viscosity. A lower bound is found using energy methods reminiscent of the 
work of Keller et al. An upper bound follows by obtaining an explicit estimate for the tractions, the 
normal component of the stress on the fluid boundary, in terms of the velocity on the fluid boundary. 
This estimate, in turn, is established using a boundary integral formulation for the Stokes equation. 
Our proof admits a generalization to other particle shapes and the inclusion of point forces to model 
self-propelled particles. 
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1. Introduction. The purpose of this work is to give a rigorous justification of 
Einstein's formula for the effective viscosity of a dilute suspension of rigid spherical 
particles in an ambient fluid. Our approach differs from those previously employed in 
the literature and is well suited to study the dilute limit for suspensions in a finite-size 
container 

Einstein [11] calculated the effective viscosity ?} of a dilute suspension of spheres 
to be 



where rj is the viscosity of the ambient fiuid and is the volume fraction of spheres 
in the suspension. By a dilute suspension, we mean a suspension where the inter- 
particle distance d is much larger than the particle size a. The formula above should 
be interpreted as a first-order calculation in an infinite asymptotic expansion of the 
true effective viscosity in powers of the volume fraction. 

Einstein's result was extended to spheroids by Jeffery in [18], who showed that 
the effective viscosity depends on the orientations of the spheroids. However, since 
the distribution of orientations is not unique in the dilute limit, Jeffery obtained 
bounds, but not an asymptotic formula, for the effective viscosity. His work was 
completed by Hinch and Leal in [TBlf2D] who computed an ensemble-averaged effective 
viscosity by including the effects of rotational diffusion, which lead to a unique steady- 
state orientation distribution. A general theory for dilute particle suspensions was 
later set out by Batchelor in [3 and Brenner in [10]. In [4l[5], Batchelor and Green 
calculated the 0{(jP) correction to Einstein's result, by including the effects of pairwise 
interactions. In [1], Almog and Brenner produced a homogenized viscosity field ?7(x) 
that includes boundary effects and agrees with Einstein's result when the volume 
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fraction is uniform. In |13fll5j . the results for disks and spheroids were extended 
to dilute suspensions of active particles (particles that propel themselves, such as 
bacteria). 

While most of the results in [DIMllinillllllSHiniliaiSQ] are based on formal 
asymptotic analysis, a more rigorous justification of equation (|l.ip is given in [19], 
using variational principles. In that work it is proved that for a suspension contained 
in a domain Q, 



as (/) — > 0, where 4> is the volume fraction. They obtain this result by placing the 
particles in a finite-size container and letting the container expand in order to avoid 
the renormalization required in |4j when considering a suspension filling M.^. The 
question of whether the same asymptotic expansion is valid in a container kept at 
finite volume when the size of the particles vanishes was not addressed. Furthermore, 
it is not immediate that the method used in [19] can be applied to the case of difi^erent 
particle shapes and self-propelled particles, as it requires knowledge of exact solutions 
to the Stokes equation for a particle inside a concentric sphere. 

Concurrently, in [5], a proof of Einstein's result including a 0{(jP) correction 
is given in the case of periodic boundary conditions. The formula is obtained by 
considering a periodic suspension of a spherical fluid drops with viscosity 77* in a fluid 
with viscosity rj in the limit 77* — )■ 00. In this paper, we consider a suspension in a 
flnite-size container and derive the formula from a PDE model with rigid particles. 

In this paper, we improve upon (jl.2l) by showing that the next correction to 
Einstein's formula for a suspension in a flnite-size container is at least of order 0"^/^. 
More precisely, we will show that there exists a positive constant C independent of 
the volume fraction <j> and the number of particles TV such that for all < ^, 



Unlike in previous work, we establish (|1.3|) assuming that the size of the container 
remains finite in the dilute limit. We model the suspended particles as spherical 
inclusions in O arranged in a regular array away from the walls of the container, 
though this arrangement is not crucial for our proof. The condition on (/) in (II. 3p 
above ensures that the spherical inclusions do not overlap. 

The method we use differs from that in |19| and admits a generalization to different 
particle shapes and the inclusion of particle self-propulsion, as modeled in [T5HT5] . We 
leave for further investigation to include collisions between particles in our model. To 
derive a lower bound for the effective viscosity, we employ energy methods reminiscent 
of the work of Keller et al [19] . To establish an upper bound, we utilize standard 
regularity estimates for solutions of the Stokes system, but we derive explicit values 
for all constants involved in these estimates. A key step in the proof in fact consists in 
bounding explicitly the tractions at the fluid boundary in terms of the velocity on the 
fluid boundary, which in turn we achieve via boundary integral equations in terms of 
the Green's function for the Stokes problem. These bounds have the same behavior 
in powers of (f>, namely (j)^^'^, thus implying (|1.3[) . 

We should remark that the cases of moderately {d « a) and densely {d <^ a) 
packed suspensions are not as well understood, though some rigorous mathematical 
results exist in both cases. Moderately-packed suspensions are studied in [53], where 
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numerical results for spheres are presented, and in |21) . where arbitrarily-shaped 
particles are studied using two-scale methods. Densely-packed two-dimensional sus- 
pensions of disks are studied using network approximations in J7^. 

The paper is organized as follows. In Section [O] we introduce the PDE model 
for a suspension of spheres and define the effective viscosity. Next, in Section 11.21 
we outline the procedure for proving equation (jl.Sp . Section [2] contains the proof 
of a technical lemma on the continuity of solutions to a boundary integral equation 
related to Stokes equation. Finally, in Section [31 the lemma is applied to prove the 
main result of the paper. Theorem 13. II 

We conclude the Introduction with some notational conventions, fl denotes a 
generic (bounded) domain in M'^ with smooth boundary. Throughout, we write || • 
llo II ■ Wl^q) (or || • \\(L\n))3 or || • || (L2(n))3x3 , as appropriate). -ff"(f2), a G 
M+, denotes the standard L^-based Sobolev space of functions on a smooth bounded 
domain, while the space (d^) of functions on the boundary is obtained via the usual 
trace relations. In addition, we define C^(f2) := {/ € C°°{n)\f{x) = for x G dVl}. 
If (T is a 2 tensor and n a vector, crn is their contraction: (crn)* = crj , while : 

denotes the standard inner product between matrices, that is: A : B := Tlij ^] ^l- 
Finally, we write / for the identity matrix. 

1.1. Preliminaries. Let 17 C be a bounded domain with a C°° boundary 
dil. Wc will model a suspension of N neutrally buoyant, passive spheres of radius 
a centered at in an ambient fluid of viscosity r] using the Stokes equation along 
with balance of forces and torques at the boundary of the spheres: 



The use of the Stokes approximation to the full Navier-Stokes fluid equations is 
justified by the small scales and low velocities (or equivalently small kinetic en- 
ergy) inherent to the problem, so that the Reynolds number is typically very small. 
Above, u is the fluid velocity, p is the pressure, Vlp '■= ^\ ^iB^ is the fluid do- 
main, g M.^ and G M'^ are the rigid translational and rotational velocities of 
the Zth sphere respectively, a := —pi + 2rje is the stress (here / is the identity ma- 
trix) and e :~ ^ (Vu + Wu^^ is the strain or symmetric deformation gradient. On 
the boundary of the domain 917, we impose a linear profile for the fiow u = ex, 
where e is a constant, symmetric, and traceless 3x3 matrix. The existence of so- 
lutions {u,p) G {C°° {p.F)Y is a standard result (see, for example, Theorems IV. 1.1 
and IV. 5. 2 in [T^). As in [TT], we will study the case when a ^ d, where d is the 
minimum distance between spheres. To simplify some calculations, we will consider 
the case where H. = B{0,R) and the spheres are arranged such that their centers 
lie on the points of the array Z'^ n B{0, R — 1) and, in particular, then d — 1. 

Remark 1.1. This geometric set up is not crucial to our method. In particular, 
our estimates are valid for arbitrary N and d as a ^ 0, so that we can assume that 

(X N d^ is fixed throughout. 

We next give a definition of effective viscosity in terms of the kinetic energy in 
the system. We denote with 
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the rate of energy dissipation in flp. Let u'* — ex. Then, we define the effective 
viscosity i) as the uniform viscosity in 57 such that the hnear flow u'' = ex dissipates 
as much kinetic energy E{u^) :— 2fi\^l\e : e as the flow u in the suspension, that is: 



1 

2'\n\e 



(1.6) 



where \ is the volume of 57. Therefore, f) is a function of the fluid flow w, ahhough 
we do not exphcitly show this dependence. 

1.2. Outline. We briefly outhne the key steps in proving our main result, the 



validity of estimate (|1.3p for 



4Tr 
3 



(a/d)^ < I (see Theorem [??T|) . This condition 



on the volume fraction ensures that the particles in suspension, modeled as spherical 
inclusions, do not overlap (recall d= 1). 

It was first shown by Einstein in [Tl] that the first-order correction ^rjcj) to the 
ambient viscosity can be formally obtained by replacing u in equation ()1.6|) with a 
dilute approximation u'^ (which will be defined explicitly in equation p.24|) ) and a 
with the corresponding stress a'^. The dilute approximation u"^ is given by superim- 
posing the Stokes solutions in the case of a single inclusion, one for each particle in 
the suspension. Interactions between particles are completely neglected in this ap- 
proximation. Our first step is to show in Lemma [33] that, using the dilute solution u'^ 
instead of u in equation ()1.6|) leads to the following bound on the effective viscosity 
t): 



1 



(i^ <fi + Cr^cjf'l^ 



for a positive constant C independent of 4> and N . Next, in Lemma 13.61 we prove an 
upper bound with the same order in (p by calculating the effective viscosity for 
(denoted 7}(m+)), where solves 



Vj3+, V 



ex 
ex 



X <^^lF 
X e dB^ 
X e dn. 



Even if we do not have an explicit form for li^ at our disposal, in Lemma 13.71 we 
obtain the following estimate on fj{u^): 



f}{u^) — ?7 1 + 



< 



1 



2e : e\n\ 



a^n{x - x^)dS 



as' 



5/3 



by utilizing the partial differential equation (PDE for short) that := — u'^ 
satisfies (see Equation p.29p ). Above, is the stress corresponding to the fiow u^, 
n the unit outer normal to the boundary, and a^n are the corresponding tractions. 
Since oc N , we have 



f)(M+) - ?7 1 + 



The last step consists in estimating the tractions, which we do not know explicitly. 
To this end, in Theorem 12. II we show that 



\a^n\ 



dn, 



<C||u^| 
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for C a positive constant depending on the domain Q,p. The most laborious part of the 
paper concerns determining an explicit upper bound for C . which entails bounding 
the pressure at the surface of the spheres in terms of the velocity. The procedure 
for doing so is outlined in Section [2] and the actual bound is given in Corollary 12.31 
Finally, ||"^||^i(af2^-) can be calculated from the PDE it satisfies. This estimate and 
the proof of the validity of equation (|1.3p (Theorem 13.11) are given in Section [3l 

2. Estimating tractions on the boundary. As outlined above in Section [TT^ 
to prove (jl.3l) we shall need to estimate the tractions t = crn on the fluid boundary 
d^p in terms of the boundary velocity values. These can be related via the following 
boundary integral equation: for ^ G dU,p (see, for example, [26 ), 



\<0 



ant 



(2.1) 



where Q — {Sm]) ^ [l /r + xx/r^) , is the Oseen tensor, V = {in) ^x/r^ is the corre- 
spending pressure, and Yiijk ~ ^VjSik + V [ 'gi^ + "57^ ) ^'^ corresponding stress. 
The solvability of (|2.ip for the tractions in Sobolev spaces is known. 

Theorem 2.1. Let f2 C M'^ be a hounded domain with a smooth boundary T = 
dfl. Let u he the solution of the Stokes equation in VL with Dirichlet boundary data 
g G (-ff"(r))"^, OL> ^, satisfying the compatibility condition g ■ ndx = 0. Then , 
there exist tractions t := an e {H"^^{T))^ such that 

l|i||(ff=-i(r))3 < C'll3ll(//°(r))3, 

where C > is a constant depending on 51, but independent of g and t. 
We omit the proof of this Theorem. 

In the remainder of this section, we will be concerned with finding an explicit 
upper bound for the constant C in Theorem 12. II above for the case a — 1. 

To this end, we will combine four inequalities, for which we evaluate constants ex- 
plicitly. First, in Section I^TT] extending the divergence- free condition to the boundary 
for regular u and a density argument gives 



\\t\\dn. 



< 



du 
dn 



du 



(2.2) 



where r', i = 1,2, are orthogonal unit tangent vectors to the surface 80, p. Next, we 
estimate the pressure and the normal derivative appearing on the right side of 
In Section [2?2] we then show that 



Ci 
where 



dn 



< 



i=l 



du 



C C 
+ — Ibllaoi. + Ci\\u\\dnj, H Ibll^F, 



(2.3) 



Ci = V2/2 - V(rf/2-a)-Ml + (4<52)-i)<54, 



^3 - + ^J{d/2 - a)-i (1 + (4<52)-i) <53 



(2.4) 
(2.5) 
(2.6) 
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C5 - V^2/(d/2 - a), 



(2.7) 
(2.8) 



for all di, 62, S3, > 0. This estimate follows from standard integration by parts 
once an appropriate extension of the normal vector n to dflp to a function on flp is 
constructed. In Section [2731 we next show that there exists C > such that 



\\P + C\\n, < {l + V2)P\\p + C\\ 



where 



1/4 



P =2i/Vl + 6/d (d/2 + fl-i + 72) 

+ 2-i/Vl + 3(i?-d/2)-i(l + i?/3)'/* . 



(2.9) 



(2.10) 



This estimate is obtained construction an explicit extension p of the boundary trace of 
p on dH. F to a neighborhood of each 95' and dH. to obtain a satisfactory dependence 
oi /3 on rip- Finally, wc show in Section 12.41 that 

du 



377, 



3?7 
R ' 



u\\dB{Q.,R) +V 



dn 



(2.11) 



This estimate is derived from a boundary integral equation for the pressure on the 
fluid boundary. The norm of the associated integral operator can be estimated via 
the Fourier transform by stereographically projecting the sphere onto a plane. Some 
of the calculations are simplified utilizing the quaternion formalism and the Clifford 
algebra ^3. 

Combining inequalities ()2.2|) . ()2.3p . (|2.9p . and ()2.11|) . we obtain (we recall that p 
is defined up to an arbitrary additive constant): 

Theorem 2.2. Let < 2a < d < R and let u,p solve 

TjAu = Vp, V-M = X e^p 

g e H^r\L°°{dnF) xedQp, 

where Jg^ g ■ ndS = 0, 17f := B{0,R) X^b, = [ji(^AB{x\a), e dl^ , and 



A = {/ G Z3 



I — a < i?} . Then there exists a normalization of p such that 



10 



6 , C4 
l_ 3 

R 2V2 
10 



II II 60 1.. ., 

WuWan^ +'7y-||w|lu,aB' 



- + 25V2 
4 



2 

E 

1=1 



du 



(2.12) 



where C4 is given in equation (j2.7p and 5k are chosen as in equation (|2.16p . 

Proof. Combining inequalities p. 21) . (12.31) . and (j2.9l) and renormalizing p such 

that 



llPl 



In^ < (1 + V2)/3||p||an, 



yields 



90j 



< 



(^3 + (1 + V2)C5/3 

2 



■»7^||u||ao^,- 



(2.13) 



(2.14) 
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Combining inequalities (|2.11l) . (|2.3p . and (|2.13p produces 
Man, [l - Cfi (C3 + (1 + V2)C5f3)' 



<^\\t\\an, +v 



3_ C4 
Ci 



3 C . , 

\\u\\anp +??-||u||u,9B' 



9u 



(2.15) 



dnf 



We wish to use inequality ()2.15p to remove ||p||aojj from equation ()2.14p . In order to 
do this, we must choose Sk (in Ci) such that 



Cfi = 



1 - (Cs + (1 + V2)C5(3)] \l + Cf' fCa + (1 + 72)^5/3 



< 2, 



which, in turn, requires that Cf^ (C3 + (1 + \/2)C5/3) < 1/3. This can be achieved 
by choosing 



<5i 



-2a 



2^ (1 + ^2)2/32 ' 



1 i?2(d-2a) r _ r,8 r 

?ni 452+1 ' ''4 - ^ O3. 



(2.16) 



Then Ci = 2V2, C2 = 32 + 1/^2, Ce - 3/5, and r/^ (C3 + (1 + 72)^5/3) = i. 

Combining equations p.l4p - p.l6p . we get the desired result. 

□ 

To apply this to the error estimate, we will use 

Corollary 2.3. Letu,p solve equation (jl.4p . Assume, without loss of generality, 
that d = 1 (and hence R cx N^^^). Then 3C > Q, independent of a and N such that 



iWdn, iV^/^ 1 

-^r- <^^Man, + -||"||u,9B' + 2^ 

2 — 1 



du 



an. 



2.1. Estimate of • In this section we derive the first of four inequalities, 

estimate (|2.2p on the tractions at the boundary oi ilp- In what follows, we fix an 
orthonormal frame {Ti,T2,n} on dQf, where n is the unit outer normal and are 
tangent to the boundary. 

Lemma 2.4. Let u,p solve 

( T]Au = Vp, V • u = X G^If (9 t7\ 

\u = geH\dQF) xeOnp, ^ ' 

where J^^^ g-ndS = 0. Then 




(2.18) 



Proof. The Lemma is an immediate consequence of the divergence-free conditions 
for smooth solutions to (12.171) and follows by density for u S {H^^^{flF))^- Let 
gm € (C°°(9r2i?))^ be a sequence of functions satisfying Jq^f 9"^ ■ ndS = that 
converges to g in {H^{dilF)Y ■ Let Um,Pm solve equation (|2.17l) with data gm- By 
Theorem IV.5.2 in [12 , e {C°^{nF))^, and therefore V • m,„ = on OQf by 
continuity. Consequently, 




(2.19) 
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Then, regularity results for the Stokes equation (see Theorem IV. 6.1 in [T^]) and the 
trace theorem give 

||Vu - VUmWdilF < Ci\\u - Um||(_H-3/2(n,.))3 < ^211^ - gm\\ {H^dnF))^ ■ 

That is, both the normal and tangential derivatives of Um at the boundary converge 
in to those of u. We can hence conclude by passing to the limit m — >■ oo in (|2.19p . 



Theorem 2.5. Let {u,p) solve equation ()2.17p . Then, 



du 



dn 



du 



dri 



Proof. Applying Lemma |2.4[ we have 



\t\\ <||pri|| + r/ll Vw • n + Vu • n\\ < \\p\\ 



<\\P\\ 



du 



dn 







du 




du 


p\\ 




dn 


+ V 


dn 



du 



dr^ 



Edu 



2.2. Estimate of 



du 1 1 
dn WaSlf 
du I 



In this section, we tackle the second inequality 
2.3p . In order to bound Hf^H^Q ' following result. 

Theorem 2.6. Let < 2a < d < R € R and B{0,R) \ VIb, where 

= \Ji^^B{x\a), e dZ^, and A ^ {l e I?\ |a;'| + | - a < R] . Let {u,p) solve 
fUJn with g e H^dVLp). Then, for all 6 := ((5i, (52, ^3, > 0, 



du 
dn 



< 



avti. 



1=1 



du 
97* 



an, 



— Wphi-iF + Ci Msuf + —WpW: (2.20) 

77 T] 



where Ci{5)~Cc,[5) are given in equations 

Proof. We again smooth out the data g first and use density arguments. Let 
g-m g (coo^^^^))3 be such that /gjj^ • ndS* = and ^ 5 in (H'^{dnF)f. Let 
{um,Pm) be the solution to (|2.17p with boundary data gm- Then, m™ G {C^ (Qf))^ , 
pm g c°°(nF), and u in {H^/^{nF)f, ^ p in H^I'^iVLF) (see, for example. 

Theorem IV.5.2 [H]). 

We next construct an extension of n e {C{dVlF)Y' to a function M e {C{VIf))'^ n 
{L'^ip.F))^ by extending it radially in the vicinity of all spherical inclusions and the 
spherical boundary. Specifically, let M € {C{VlF)Yf^{L'^{^F)Y be such that A/" = n on 
d^F- Let A/" = on dB (a;', f ) for ^ G A and extend TV inside B {x\ ^)\B{x\a) such 
that each component is linear on any line emanating from . Similarly, let A/" = on 
dB (0, i? — I) and extend Af inside B{Q, R)\B (0,R— |) such that each component 
is linear on any line emanating from the origin. Where it has not yet been defined, 
set Af = 0. It is easy to check that Af has weak derivatives in L'^{^lF)<^L°°{^lF), that 



IIA^II 



(L~(Of)) 



3 < 1, 



|VAA|| 



< 



d/2-a ' 



(2.21) 



and V • A/" > almost everywhere. 
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By the definition of A/", we can write 



dn 



dS = 



dx 



du" 
dn 



■dS = 



dx 



du" 
dn 



■dS 



- [ V • [iVu"W) Vu™] dx. 



Performing the differentiation, using the Stokes equation, and integrating by parts, 
we get that 





du"' 




/ 


dS^ f 


Janp 


dn 


Jn 



If Vu"" ■.Vu'^dS -I [ (V • AA) Vu" : Vw™ da; 
+ [ Vu"V7V : Vu'" dx + - [ (Vm"7V) • Vp" dx. 



Integrating by parts once more and using the fact that V • = 0, we get that 

2 





du"' 


/ 


dn 



dS : 



1 
2 

+ 



Vu"'VJ\f : Vu"' - -Vu" : (WAff p" 



dx 



- f Vu™ : nnp'^dS -]- j {V ■ N) Vm™ : Vu™ dx 



=:/i+/2+/3+/4. (2.22) 

Since V • A/" > a.e., we can neglect I4. We can decompose Vu™ = Wu"^nn + 
J2i=i ^u"^tW^ to estimate 

1 



du^ 



dn 



1=1 



di 



dri 



dS. 



(2.23) 



We can use Cauchy's inequality to estimate, for all 5i > 0, 

2 



h < 



1 

4Sl 



dQ.F 



du" 



dn 



n 



Applying Lemma [2^ we get that 





du"' 


1 


dr^ 



dS+^\\pn\lnF 



dS+^\\pn\laF 



(2.24) 



We can also use equation (|2.2ip and Cauchy's inequality to get, for all 82 > 0, 



I2 < 



1 



<- 



d/2 - a 
1 

'd/2 -a 



E 

3,k 



du"' du''^ 



4^2 



dxj dxk 



dx 



2J 



du 



dx 



dx 



(2.25) 
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Integrating by parts and using the Stokes equation, 

J<rtp Jdn^ on 

= -- [ u"'-Vp"'dx+ [ vJ^Vu'^-ndS. 

V JQ.F JdQ.F 

Integrating by parts once more and using the fact that V • — 0, we get that 



n\L^- u^-np"^dS+ ^'"-V 



Using Cauchy's inequality twice more, we get that 



dn 



(2.26) 



Combining estimates (|2?22|) . (|2?23l) . ([Z^M]) . H^^ . and (fZ^Mj) . and taking the hmit as 
TO — >■ oo (all terms converge to the desired quantities due to Stokes regularity-see, for 
example. Theorem IV. 6.1 in [T2]-and the trace theorem), we get equation (|2.20p . 
□ 



2.3. Estimate of We now turn to proving the pressure estimate 

We will normalize p, which is defined only up to an additive constant, in order to 
optimize constants in that inequality. 

Since best constants in trace and Sobolev inequalities are difficult to obtain in 
general, we proceed as follows. We first construct an explicit extension p £ H^/'^iTlF) 
of the trace of p in Lp'idVlF), for which we have good control on the trace constants. 
Then, since p is harmonic from Stokes equation, p = p po, where po solves 

Apo = Ap X € f^F fr, r,7N 

Hence, < ||p||Hi/2(nF) + IIPollfiF- Then, WpWuf < (1 + C'i)C2||p||aoF, where 

Ci, given in Theorem 12. 81 is the constant in ||po||aF < C'i||p||ffi/2(siF)- and C2, given 
in Theorem 12. 15[ is the constant in ||p||^i/2(Qp) < C2||p||aoF- 

Lemma 2.7. Let p G f/^/^ n L°°(17f) and let po e H^^'^iflp) be the solution of 
equation (j2.27p . Then 



Ibollo, < IIpIIo, + llPlli-on,) ■ (2.28) 



Proof. Let p„ £ C^iilp) be such that p„ — > p in H^/'^{Q,f) and let p^ be the 
(smooth) solution of (|2.27l) with datap„. Then, by standard regularity results (see, 
for example, [21]), pI ^ Po in H^/^np). Let £:={</.€ C^inF)A\M\nF = !}■ 
Using the dual characterization of the L^-norm and the density of in gives 

lb° Ibj, = sup / Apn(f>dx = sup <^ - / Vpn ■V(f>dx\ 

<pe£ JUf 06£ L JUp ) 



snp< / pnA(j)dx- I Pn-^dS 
0e£ VJuf JdnF 
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dn 



dS. (2.29) 



We next note that for any (f> in the admissible class £, one can find a sequence 4>k & £ 
such that > and — >■ |^ on dVLp in L^{dVLF)- Therefore, we can assume 
1^ > without affecting the supremum. Therefore, 

WpIW <\\Pn\\n^ + \\Pn\\L^(Qn^) sup f ^dS 
= \\Pn\\n^ + \\Pn\\L^i^Qn^) sup / A(/)da; 



Since po and p„ p respectively in L'^{VIf) and L°°(951f), estimate ()2.28p 

follows. □ 

Theorem 2.8. Let p e H^/^ n i°°(r2i;^) and po e &e the solution of 

equation (|2.27p . T/ier* f/iere exists a constant C > smc/i i/iat ||polkiF — V^Hp + 

Proof. We can add an arbitrary constant to p without changing (j2.27p , and hence 
its solution po- We hence can substitute p+ := p + 3||p||i,oo(f^^) for p in (|2.28p to 
obtain: 

IIp+IIL ^ 



> 

-2 



(p2 + 6p||p|U^(a,) + 9|!p||i„(o,)) dx > WPWl^ + 3 \nF\ ||p||i^(o,) 

]^ / s 2 I 



□ 

In order to prove Theorem I2.15[ we will need the following Lemmas. We include 
proofs for the reader's sake. 

Lemma 2.9. Let p e H^/'^{dB{0, R)) and let p'^ e H^{B{0,R)) solve 



Ap^ = X e B{0, R) 
p'^p xedB{0,R). 



..^R^n m (2-30) 



Then ||p^||B(o,i?) < y f l|p||aB(o,fl) and || Vp^||s(o^fl) < \p\m/2(^oB(o,R))- 

Proof. We first recall that for any s > 0, linear combinations of spherical har- 
monics are dense in H''{S^), where = dB{0,l) is the unit sphere, and that an 
equivalent norm is given by 

/ 00 n \ 1/2 

11/11^.(5.) = 5] 5] A^1G„.N , (2.31) 

\ri=0 fe=-n / 

where A„ — n (n + l), and Ank is the L^ projection of / onto YJI;, a spherical harmonic 
of degree n (we refer, for instance, to \22\]. By rescaling a similar expansion is valid 
for p € L^(9_B(0, R) for any fixed R> 0. Moreover, if we set p^{u}) := p{Ruj), uj e 5^, 

\\p\\HHdB(o,R)) = R^'' Wp'^'Wh^^s-)- (2.32) 
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By density, we can assume therefore that p is a finite hnear combination of spheri- 
cal harmonics: p = X]!^=o ^2=^n ^nk{R) Yn.k- The solution p^ of the Dirichlet prob- 
lem (I2.30|) can then be written (see for example [221 Theorem 2.10, p. 145]), using 
the homogeneity of the Laplace operator, as p^ = X]n=o X]fc=-n ^ " ^"fc(-R) ?'"'^n ) 
which is in fact an harmonic polynomial in the whole of R'^. 

Next, we observe that, if (r, 9, (p) are spherical coordinates, then 



N 



n—0 k— — n 



where n is the unit outer normal to dB{Q,R). Integrating by parts and by the or- 
thonormality of the spherical harmonics, we hence have: 



|V/| 



B{0,R) 



dB{0,R) 



/ N n \ N n 

\n—Oq— — n / n— g— — n 



(2.33) 



|p|_f/i/2(as(o,_R)) = 



|G;m(i?)p3^, 



\\P\\dB{0,R) 



Additionally, from (P3T|) -(P3 ^ . 

T.f=,T.L=-i\Gim{ R)\'' R\ and 

YZo EL=-/ Gj^^l{l + l)R. Since + < 1 for Z > 0, we get the de- 

sired estimate for p a finite sum of spherical harmonics. The estimate for any given 
p G H^/'^{dB{Q,R)) then follows by density, using the continuity in Sobolev spaces of 
solutions to (11301) • □ 

Lemma 2.10. Let p e H-^/'^{dB{Q,R)) and let p^ solve equation Then 

p' e L^{BiO,R)) and \\p'\\Bio,R) < \/l\\p\\H-^/HdB(o,R))- 

The proof of the Lemma is very similar to that of Lemma 12.91 

Lemma 2.11. Let p e L'^{dB{0,R)) and let p'^ e H^/'^{B{Q,R)) solve equation 

^M- Then, ||p^||ffi/2(B(o,i?,)) < {\)~' (1 + f )^ \\p\\dB(0,R)- 

Proof. The result follows by interpolating the estimates in Lemmas 12.91 and 12.101 
(see e.g. Theorem A. 1.1 in p7] ). □ 

Lemma 2.12. Let p £ L^{dB{0,a)). Then there exists an extension p^ € 



i/i/2(B'=(0, a)nB(0, d)), where d> a, such that 



P 



\dB{0,a} 



— p satisfying 



P ll//i/2(B':(0,a)nB(0,d)) 



< 



2i (d+i + V2)^|b||as(o.a). 

Proof. The proof is similar to the combined proofs of Lemmas I2.9H2.111 □ 

To obtain an appropriate extension of p in Q.f from Lemmas 12.91 - I2.12[ we 

introduce the cutoff functions 




i-(5^)^ 



, a < r < d, 
r > d, 

, d<r <R, 
r < d. 



(2.34) 
(2.35) 
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One can easily calculate that 

3 

liC^||L~(B<=(0,a)) < 1, llVC^||L~(B<=(0,a)) < 

J (2.36) 

IIC^I|L~(B(0,i?)) < 1, l|VC^||L~(B(0,fl)) < ^■ 

The next two results follow easily from (|2.36p . 

Lemma 2.13. Let V{a,R) := B=(0,a) n B{0,R) with a < d < R and G 
H^/^{V{a,R)). Then p^C^ e H^/^{B%0,a)) with 

|lp^C''l!ffi/^(Be(o,,)) < (^1 + ' WP^'WHU^iViaM))- 

Lemma 2.14. Let p^ e H^/^{B{Q,R)) with d < R. Thenp^C^ e H^/^{B{0,R)) 

1 

with \\fC^\m/2^B{o,R)) < (l + l)' \\p^\\m/2{B{a,R})- 

We are now ready to state the extension theorem for p. 

Theorem 2.15. LetO<a<d<ReR and B{0,R) \ Qb, where VIb = 

[ji^j^B{x\a), cc' e dZ^, and A ^ {I e \x^ \ + ^ - a < R} . Let p e L^idVLp). 
Then there exists an extension p S H^^^{U,f) such that — p, which satisfies 

I1pII//i/2(j7j^) < /3||p||90f' where (5 is given in equation (I2.10p . 

N 

Proof. Define p{x) — — — a;') +p^{x)C^{x), where p^ and p^ are 

1=1 

given in Lemmas 12.121 and I2.11[ respectively, and and are the cutoffs defined in 
equations (j2.34p and (|2.35l) . respectively. Thenp|g^ = p and, from Lemmas 12.111 - 
[2l4l it follows that 

jY 1 1 

WPWHU^n.) < E (l + ' (^^ + i + \\phBi.',a) 

11 1 

+ (i + ^ (^) ' (i + f ) ' MUio^m < mu.. 

Q 

Combining Theorems 12.81 and 12.151 '^e finally obtain (j2.13p . 

The results in this section can be extended to other shapes through the use of 
eigenfunctions of the Laplace operator on their surfaces. For many simple shapes, 
including ellipsoids, these are explicitly known |17) . 

2.4. Estimate of ||p||aap. In this section, we establish the last needed estimate 
(|2.1ip . We will do so by employing the boundary integral equation for p derived from 
the Stokes system, where we recall P(i^) = z^/(47r (see, for example, [26]): 

lpiO = -lp-y- I r{^-iy)-t{iy)dS,-rjVi-p.y. [ nV[^-v)-u{u)dS,. (2.37) 

To simplify notation, we will henceforth drop the p. v. from principal value integrals. 
Any integral involving V in this section should be interpreted in a principal value 
sense. 
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To explicitly estimate the norm of the kernel V in each integral, we employ the 
Fourier transform. To this end, we embed into C'' and thence into A^, the 8- 
dimensional Clifford algebra CZo,3 with basis 1,61,62,63,6162,6263,6163,616263. For 
simplicity, we will write 6y := 6^6^ and 6123 := 616263. Under this embedding, 
scalars are identified with multiples of 1 and vectors x = (xi, a;2, 3:3) e C'' with 
a;i6i + a;262 + X363. For an arbitrary x = xo + xiei + • • • X126162 + • • • G A3,, we let 

X := Xq- xiei - X2e2 ~ 0:363 - 0:12612 - X23623 - 2:13613 + 2:1236123, 
X := xo + XiCi + 0:262 + 2:363 - 2:12612 - 2:23623 - 2:13613 - 2:1236123, 
x^ 2:q - xld ~ 0:262 - 2:^63 - 0:12612 - 0:23623 - 2:^3613 + 2:^236123, 



where x* denotes the complex conjugate of Xi. 

Let ri,r2 C be spheres or planes. By l?{Tk,A'i) we denote the space of 
functions u : — > ^3 such that || Jp^ vix^v' (x)dSx\\ < 00, which is a Hilbert space 
with the inner product 

(/,.g)p^ P{x)gix)dS,. 

Since this definition is valid for all i^(Ffc,^3) functions, in which C°°{Tk,A3) are 
dense, we can extend this definition to the case where either / or g is a distribution. 
Under the linear fractional transformation $(0:) = (ax + h)(cx + d)^^ mapping F2 to 
Fi, we can write 

(/(2/),5(y))r, = {J{^,x)fmx)),J{^,x)g{^{x)))^^, (2.38) 

where J(<i?,x) (cx + d)||co: + d\\^^. 

We now assume T — and set (j){x) — (eix + 1) (x + 61) ^. Then, — is 
the standard stereographic projection from onto M^, identified with the span over 
M of 61,62. This projection allows one to define the Fourier transform for functions 
feL^iS^As) as: 

J-52(/)(z):=^J(V,z) / 6^W-)''^W>7(V^,o:)/(x)d5., 

where J{ip,x) = 4 (0; — ei) ||o: + ei||^'^ and {x,y) :— J^^^i^iV-i- -^s^ is an isometry 
from L'^{S'^,A'i) mto itself (see e.g. [28].) 

Theorem 2.16. Lett(^ {L'^{dB{0,r))f and 

TloB(o,r)t:^ f V{£.-y)-t{v)dS,. (2.39) 

JdB(0,r) 

Then, ^dB(o.r) "i^ ^ hounded operator on L^{dB{0,r)) with norm \\^dB{o,r)\\ ^ ^■ 

Proof. We remind the reader that the integral in (|2.39l) must be understood in 
a principal value sense. This integral is invariant under dilations, so it is sufficient 
to consider only the case r — 1. Moreover, it is enough to calculate the L^-norm of 
J^52 (1152 1), as the Fourier transform is an isometry. Let tm € C°° (F^ , ^3 ) be such that 
tm ^ t in L^(Ffe,^3). Then, noting that (see, for example, [5S]) V ~ ''IK^)) = 
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J{^p,^)-^Vi^- iy)Ji4>,i^)-'^, we have 



The above equality is vahd in principal value sense (i.e., pointwise outside a neigh- 
borhood of radius e about = v = — ei, since J is singular there, and then passing 
to the limit e — ^ 0). Since J — J, and ,P = J (because the components of J are 
real- valued) , this is 

= j(v^, ( j(^, ^) j(^, ^r^ {-Pim - ^(^)))^ , Ji^^, ^)tM)^, ■ 

By equation (|2.38p . we obtain 



We observe that, since J J ^ = 1, J ^(t/ijx) is bounded over the sphere and vanishes 
for X — — e^, so that J{ip, (^(^(0 ~ t^))^ is well defined as a principal value 

distribution. Therefore, this and the following expressions are well-defined. Next, we 
compute the Fourier transform of this expression: 

(n5.i„)(z) =^J{i^,z)(^J{^,Oe-^^^^-^^'^^^^\ 



Denoting convolution by *, it follows from equation (|2.38l) . 



27r 

=j(vaz)j-m2 (vio) wz))j'm^ {Jii^,mr't„rm))) (^(^))- 

(2.40) 

We now observe that, since J{(l),ip{S,))^^ — J('0,C) ioi ^ e M^, 

= - J(^,z)-lj-52 (i) (Z), 

again by (|2.38p . In addition, "P if a Fourier multiplier with bounded symbol {V{^)), 
since (see e.g. [13]), 

^R.(7'(0)(^^) = -^^, (2.42) 
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and, therefore, we can pass to the limit m ^ oo in equation (|2.40p . FinaUy, we 
combine this equation with (|2.4ip and (|2.42l) . and observe that for z G M'^, z)^^ = 
-J(V',^)|| JlV'j-z)!!"^, so that 



2 \\V[z)\\ 
^ i J{tp,z) Tpiz) J{ip,z) 

2||J(V^,z)||||V'W||||J(V,^)ll 

□ 

Theorem 2.17. Let T = [jidB{x^,a^), where dB{x"\a"') n 5B(x",a") = 
Vm ^ n. Let t G L^{dT,A3). Then, as an operator from L'^{T,A3) — L^(r,^3), 

linrll < i 

Proof. We introduce the mapping T^{f){z) := Y.i^dB{xKa'^){f){z)XdB[x^ ,a^), 
where Xf is the characteristic function of the set F. We observe that -^^^(r) = 
®/L^(9i?(a;', a') and by construction J-^ [f) — ®iJ'aB(x',a')- Consequently, J^j? is an 
isometry on L^{T,A3) as J^dB{x^,a') is an isometry on (L^(9i?(x', a')))^. Theorem 
UTithen imphes that HHrll < 5. □ 

We now turn to study the second integral on the right hand side of equation 
(|2.37p . which we rewrite making the change of variables ~ — v, as 



• / n{iy)ri^ - y) ■ u{v) dS, = • / - C)7'(C) • u(C - C) dS^, (2.43) 

where £, — dil.F :={C = C~^j ^£ d^lp}- Recall that this integral must be interpreted 
in a principal value sense. We temporarily assume u € C'°°{dQF) and note that in the 
sense of distributions Vllu = 11 Vu. Consequently, from the definition of distributional 
derivative: 



[Vc • n(c - OViO ■ u{i - c) + n(e - C) • (Ve«(e - Ovm ds^ 

[Vc • - ono ■ u(e - c) + «(e - • i^M^ - ovio)] dsc 

[V, • n{i^)V{( - v) ■ u(y) + n{v) ■ ,u(y)V{i - v))\ dS,, 

where the last equality follows by setting v ~ — Q. Since is a Fourier multiplier 
in L^, this last expression is valid by density for any u € H^^dflp)- Furthermore, 
for X e dB{x\a), W ■ n — - (or — - if n is the inward-facing normal). Thus, for 
np -.^ B{0,R)\uf^^B{x\af, 

f ^33 du 

• / nV{^ - ly) ■ u{v)dSy = - V -XIqb(x'- .a)U + -s^aB(Q.R)U + Hao^ — . 
J drip an on 

Under the hypotheses of Theorem 12.21 Theorem 12.171 then gives 



V^- / nV{S,-v)-u{v)dS^ 

dflp 



du 



On 



;2.44) 
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Finally, combining equations (j2.37p . (|2.44l) . and Theorem 12 .171 we obtain 



dn 



V 



du 



dn 



The results in this section can be extended to other particle shapes by using 
mappings from their surfaces to S^. 

3. Error estimate. With estimate (I2.12[) on the boundary tractions in place, 
we can now turn to establishing the needed error bounds for the effective viscosity. 
As in Section ll.il we model a suspension of rigid spheres of radius a, centered on the 
lattice n B{0, R — 1), that is, d = 1. The limit 4> '■= ^cl^ then corresponds 
to a — >■ for a fixed, but arbitrary, number of spheres N. Recall that the effective 
viscosity is defined by 



fi = {2\n\e:e)-^E{u), 



(3.1) 



where E{u) :— 2rj J^^ e : edx, e — ^ (Vu + Vu"^), and e is the constant, symmetric. 



trace free matrix describing the boundary conditions in equation ()1.4p . We will bound 
77 below and above by quantities which have the same asymptotic behavior to order 



Theorem 3.1. 
such that yO < 6 < 



There exists C > independent of a, N , and hence of <j) '■- 
^ the effective viscosity f) satisfies 

5 



3 " ' 



?7 — 77 1 + 



< C7;03/2 



Remark 3.2. The upper hound for (f> comes from the condition that the spherical 
inclusions do not overlap. 

Proof. The proof of Theorem 13.11 is an immediate consequence of Lemmas 13.51 
3.61 and 13.91 given below. □ 



3.1. Lower bound for the effective viscosity. In order to obtain a lower 
bound for 77, we will use a similar technique to that used in [TH]. As in [T^, we will 
bound E(u) from below, but we will use the dilute approximation to u (the solution 
of equation (jl.4p ) given by 



N 



'{x) := ex + ^[u^{x - x'^) - ex] , (3.2) 



1=1 



where solves 



7/Awi = \/p\ V • = a; e ]R3 \ B{0, a) 
u^=0 xedB{0,a) (3.3) 

^ ex a; — >■ 00 



and is given by — ex + u^' , with 



'2 

The corresponding pressure is given by 



u '{x) = —ex- — -J- — x{e : xx)- - — — - — ) . (3.4) 



pi(x)=pi'(a;) = -5aV:^- (3-5) 



18 



B. M. HAINES AND A. L. MAZZUCATO 



In order to obtain error estimates independent of the number of inclusions N, we 
will need to regularize various sums. We begin with two lemmas. 
Lemma 3.3. For all x e B ^0, and allQ< p <oo, 



1 



E 

zGZ3\{o}ns(o,p) 



Lemma 3.4. For all x e B{0, the 



1 



< (876 + 504y3)7r. 



(3.6) 



E 

zez3\{o} 



< 



27r (^159 + 92V3 + log 4 - log(7 - 4^3) 



(3.7) 



Both estimates follow by replacing the sums with integrals and using spherical 
coordinates to bound the integrals. The convergence of the integrals was checked by 
hand and the explicit numerical bounds were computed by Mathematica. 

Next, we show below that we can use u~ , defined in equation (13.21) . to obtain a 
lower bound for 77. 

Lemma 3.5. For fj defined in equation (j3.ip . there exists C > independent of 
(j) and N such that V0 < |, V + ^(f) < V + Crjcj)^/^. 
Proof. The proof is a calculation. First, note that 



E{u - U-) = Eiu) + E{u~) - At] 



e dx. 



(3.8) 



Now, integrating by parts and taking advantage of the boundary conditions on u, we 
have 



277 



; dx — 



a : Vm dx = 



N 

E 

1=1 



/ an 

JdQ 

(v'+w' X {x-x^))dS = 



dS 



a n ■ eX' 



dS, 



on 



(3.9) 



since the particles are force and torque-free (that is, /^^i cr n dS = and J^^i a 
{x-x^)dS = 0). Now, 

E{u^) — 2 / n ■ exdS = / : dx — 2 / a^n ■ exdS 

Jon Juf. Jan 

-'ds+y I 



a n ■ exdS 



I en 

--■.h+h+h. 



a n ■ u 



a n ■ u dS 



(3.10) 



First, we can write a — 2r]e + a ' and integrate 1\ by parts to obtain 

l\ — — \ a^' : edx + / a^'n ■ exdS — 2?/ / en ■ exdS 
Jhf , J SB' Jan 



27/e 



[ e-'dx + yf a-'n-exdS 

Jqf , JOB' 



Ion 
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Jdn I JdB' 

+ / a^'n • exdS 
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e : e. 



Some of these integrals above can be evaluated explicitly (see e.g. [5]) to give 
/,^_^,:efl + ^0)-2,y / e:u-'ndS 

+ ^^e: f [a^'{x-x"')nx^2rju^'{x~x"')n\dS 



(3.11) 



dB' 



Now, by integrating by parts and recalling that n is the inward- facing normal to dB^ , 
we have 

f ai'(a; - = - V e : / a^' {x - x"') dx ^ 

m^l J^B^ .n^l 

- 277 Ve: / e^' {x - x"") dx ^ 27]^ e : [ u^' [x - x'^)n dS . (3.12) 



Combining equations (jS.lip and p.l2p gives 



MRS o _ 

ii = e : e 1 H — (p I — 2?7e : / u ndb. 

3 V 2 / ./gn 



(3.13) 



We can split I2 into two integrals by writing again a = 2rje + a ': 

I2 ^ <J^ : u^'ndS ^ 2r]e : I u~'ndS + / a^' : u^'ndS. (3-14) 
Jdn Jdn Jdn 

The first term above will cancel the second term in equation p. 131) . Now, we can 
estimate 





: u 'n dS\ 


(3.15) 


Jdn 







j ya^'{x-x^):yu^\x~x"')ndS 
Jdn , ™ 



which, using the formulae for u^' and p^', given in equations p.4p and p.Sp . respec- 
tively, and the fact that J^^^ a^'ndS — 0, is 



<c\dn\ 



a'^ inf 



ya^'ix^x^) 

E 

-nB{ 

ya^'ix-x^) 



(L~(ao))3 
1 



ez3\{o}nB(o,2Afi/3) 1^ ^1 



(L=°(an))3 



C 



0) 



20 



B. M. HAINES AND A. L. MAZZUCATO 



1 1 



E 

zGZ3\{0}nB(0,2Afi/3) 



Applying Lemma l373l and recalling that cx N, this is 



0) 



<CiV2/3 



n3 



2eZ3\{o} 



< Ca^N^/^loe-N. 



We estimate I3 in a similar fashion. First, we note that u^{x — x') = u (x) 
12m^i — x™) is constant on dBK Hence, 

Is^S^ f cr- -.u-ndS [ a- -.y^ u^'{x - x"')ndS 

I J9B^ I J^B^ m^l 

To regularize the sum above, we set for a; € Z^, 

^„ r -pi'(x)/ + r/ (Vw'^'(x) + (Vm'''(x))^) a; e Z3 \ {0} 

|_ a; = 0. 

We divide the integral further by writing = 2rjE + a^': 

/g^V/ a- -.J" u^'{x- x"')ndS 
I -^93' 

--27jS2E:[ J" u^'{x-x"')ndS + y" [ a'' : u^' {x - x"')ndS 
, Job' z^, , Jos' 

(3.17) 



(3.16) 



Given Lemma [3.31 and the formulae for u^' and p^' , given in (|3.4I) and p.Sp . respec- 
tively, we can estimate in a manner similar to equation p.isp . taking advantage 
of the fact that Jg^, ndS ^ 0: 



27]S2E: [ S2u^'{x-x"')ndS 



<C?7a^y las'l inf 



E 



c 



2eZ3\{o}nB(o,2Ari/3) 



E 

2ez3\{o}nB(o,2Afi/3) 



1 1 



L°°(aB(0,a)) 



L°°(aS(0,a)) 



(3.18) 



<Cr]Na^. 

To estimate /|, note that, by the Cauchy-Schwarz inequality and the fact that J^^, cr^'r 
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0, 



l7"M < 

r3 — 



V / cr-' : V u^'ix - x"')n dS 

I JSB^ m^l 



<Ca^ llCT^'nll „„, iV^/^ jj^f 



E 



c 



zez3\{o}ns(o,2Ari/3) 



E 



zeZ3\{o}ns(o,2Ari/3) 



X — zr z 



aB(0,a) 



9B(0,a) 



IU,C 



Proceeding as in equation p.lSp . since Jg^, a^' (x — x^)ndS = 



^a^\x-x^)^ 



< 



UidB' 



J2[a"{x-x'')-a^'ix'')] 



Making use of the formulae for u^' and p^' , given in equations p.4p and p.Sp . respec- 
tively, and Lemma [3. 4[ we get that this is 



y{x~x^)~a'-'{x^) 



(L°°(u,aB'))' 



E 



1 



2ez3\{o} 

The bounds ((3T7H3TT91) imply that 



\x — z\ 



< Cr]N^^^a\ (3.19) 



L°°(aB(0,a)) 



I/3I < Cr]Na^ 



(3.20) 



We recall that \n\ = A/SttR^ cx N. Then, p.ip and previous estimates combined 
finally yield: 

for some constant C > dependent on u, but independent of N and (f>. □ 

3.2. Upper bound for the effective viscosity. We will next show that u+, 
solution of the boundary value problem 



77AM+ = Vp+, V • M+ = xe^F 

X G dB^ 
X e dn, 



u+ = ex'' 



(3.21) 



provides an upper bound for the effective viscosity. 

Lemma 3.6. Let fj he defined in equation p. 61) and let 



fi{u+) -.^ {2\n\t : e)~^E{u+). 



(3.22) 
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Then, f) < fi{u^). 

Proof. We begin by writing 

E(u+ - m) = E{u+) + E{u) -2r] e+ : e dx, 

Next we observe that, since the particles are force and torque-free, 



E{u) = 277 



e : e 
an ■ ex 



dx = 



I dx = 



a : Vu dx 



/ an ■ ex dS + y~] / an ■ (i;' + x (a; - a;')) dS 
Jon , Job' 



an ■ ex dS 



an 



IdB 

: a dx — 21] 



e dx. 



Therefore, E{u+-u) = E{u+)-E{u) > 0, that is, E{u) < E{u+), which by definition 
imphes i) < fi{u'^). □ 

In hght of Lemmas 13.51 and 13. 6[ it remains to prove that there exists a C > 
independent of cj), N such that for all < ^, 



r]{u+) - 77 1 + - 



(3.23) 



To do so, we introduce a dilute approximation for u 



+ . 



:=Ex + Y, - x^) - Ex] - C"^, 



(3.24) 



where solves equation (|3.3p . The constant C^, which will be determined later, is 
added to regularize the sum-that is, to ensure u"^ remains finite as A'^ — >■ 00. 
We begin with a preliminary lemma. 

Lemma 3.7. There exists C > independent of a and N such that 



ri{u ) — 77 1 + 



< 



2e : e\n 



a'^n(x^x']dS 



+ (3.25) 



where a^ is the stress corresponding to the flow := — u'^ , where u"*" is the 
solution to equation (j3.2ip and u'^ is defined in equation (I3.24p 

Proof. We integrate the right-hand side of equation (I3.22p by parts, utilizing 



2me:eJ^ 



a^ : dx 



a-nxdS + Y. ---'dS 



dB' 



where by further integration by parts we can write (recall a^ = —p^I + 2rje^): 
e:f a-nxdS^2r,e: f e+ - V . : / a^nxdS 

=2,e: / exndS + e:y:f {2,enx^ - a- nx) dS 
Jan , JdB' 



(3.26) 
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(3.27) 



We then have, reversing the orientation so that n is the outward-facing normal to 

fi{u+)=ii+ —-^ y^^ - I (T^n{x- x^)dS. 

Substituting cr+ = 2rie + cr'^' + in the above equation, and integrating exphcitly 
some terms (see, for example, [6 ) gives 



x')dS 



+ ^^—ye: I ya"ix~x^) 



n [X — x^) dS. 



We denote the last integral above by /, and we observe that, since solves 
integrating by parts, 

/ = - e : / V (V • a^'x + a^' (x - x"')) dx 
= -2Tje: f y e^'{x-x"')dx = -2r]e: f Y u'' {x ~ x"')ndS 

Jb' JdB' 

Note that, since J^^, ndS = 0, by Lemma [3731 and formula p.4p for u^', we have 



/ S2u^'{x~-x"')ndS 
< Ca^ \dB^\ inf 



E 



c 



!ez3\{o}ns(o,2Afi/3) 



E 

zez3\{o}nB(o,2AfV3) 



(L~(aB(0,a)))3 



(L°=(9B(0,a)))3 



Therefore, recalling that |r2| oc N, 



I mf^l 



X — x'') dS 



The next Lemma is a straightforward consequence of Lemma 13.71 and Corollary 
Lemma 3.8. There exists a C > independent of 4> and N such that for all 
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< 



— ?7 1 + 



<Cr]- 



/,l/3 



/,l/2 



N 



lUidB' 



02/3 



(3.28) 



5/3 



Proof. From Lemma 13.71 we can write 



fi{u+) - ?? 1 + 



< 



< C 



5/3 



where the second inequaUty follows from the Cauchy-Schwarz inequality, noting that 
|a; — x'l = a for x £ dB^ and \Q\ cx N . Then Corollary 12.31 implies that 



m(u^) —77 1 + 



< C 



T7 C 'T- an 



7,3/2 



/,5/3 



„ / a II p II EL ' II p II a 

— I II iiu,as' ' ]\[i/2 II iisof 'TV 



■qcj)- 
a' X - 



which, upon making the substitution 1 — (^^^ , is precisely equation p.28p . 
□ 

The upper bound (|3.23p will now follow from the previous lemma once we estimate 



\H^(dnf 



For this purpose, it is useful to write a PDE that satisfies: 



X e dB\ 7^ 
X G 

a; G dVl, 



(3.29) 



where 6' = X^^i "^'(^ " '^''^ " '^'^ = ^fc '^^'(^ " '^''^ " 

We will use this PDE to prove the main result of this section. 
Lemma 3.9. There exists C > independent of a, N, and hence (p, such that, 

for all < 4> < ^, equation p. 231) holds. 

Proof. From Lemma [3.81 we have that there exists a C > independent of N,a, 

and hence A such that 



fi{u+) - ?y ( 1 + -<; 



a/2 



]\fi/2 11" iiaOf 



(3.30) 



Now, using the formula for u^' in equation (13. 4[) . selecting such that the 
infimum 



inf 



(L~(aaF))3 
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is achieved and applying Lemma 



E 



zez3\{o}ns(o,2Ari/3) 



1 

w 



L°^{aB{0,a)) 



(3.31) 



Similarly, 



(3.32) 



E 



<C7Vi/2||Vu^| 



C2N 



Combining equations p.30H3.35)) . we have the desired result. 



(3.33) 
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